Thermoballistic spin-polarized electron transport in paramagnetic semiconductors 
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Spin-polarized electron transport in diluted magnetic semiconductors (DMS) in the paramagnetic 
phase is described within the thermoballistic transport model. In this (semiclassical) model, the 
ballistic and diffusive transport mechanisms are unified in terms of a thermoballistic current in 
which electrons move ballistically across intervals enclosed between arbitrarily distributed points 
of local thermal equilibrium. The contribution of each interval to the current is governed by 
the momentum relaxation length. Spin relaxation is assumed to take place during the ballistic 
, electron motion. In paramagnetic DMS exposed to an external magnetic field, the conduction 

■ band is spin-split due to the giant Zeeman effect. In order to deal with this situation, we 
extend our previous formulation of thermoballistic spin-polarized transport so as to take into 
account an arbitrary (position-dependent) spin splitting of the conduction band. The current 
and density spin polarizations as well as the magnetoresistance are each obtained as the sum of 
an equilibrium term determined by the spin-relaxed chemical potential, and an off-equilibrium 
contribution expressed in terms of a spin transport function that is related to the splitting of the 
spin-resolved chemical potentials. The procedures for the calculation of the spin-relaxed chemical 
potential and of the spin transport function are outlined. As an illustrative example, we apply 
the thermoballistic description to spin-polarized transport in DMS/NMS/DMS heterostructures 
formed of a nonmagnetic semiconducting sample (NMS) sandwiched between two DMS layers. 
We evaluate the current spin polarization and the magnetoresistance for this case and, in the 
limit of small momentum relaxation length, find our results to agree with those of the standard 
drift-diffusion approach to electron transport. 

S ' PACS numbers: 72.20.Dp, 72.25.Dc, 75.50.Pp 

O ' I- INTRODUCTION 

The electrical injection of spin-polarized carriers from magnetic contacts into, and their subsequent transport across, 
nonmagnetic semiconductors (NMS) constitute outstanding issues in the field of semiconductor spintronicsJ^'^'^'^ 
^ While spin-based semiconductor devices still await realization, major progress has been achieved over the past years 
in the understanding of the conditions and mechanisms governing polarized-carrier injection and transport. 

■ The low injection efficiency observed for ferromagnetic metal (FM) contacts found an explanation^, in the large con- 
, ductivity mismatch between the contacts and the semiconducting sample. To circumvent this obstacle, the introduc- 
tion of spin-selective interface resistances was suggested^i^ and analyzed theoretically in some detaili^ i-'^^'^^i^^'^'^i-'^'^'^^i^^ 
A considerable enhancement of the injection efficiency due to externally applied electric fields was predictedj^^iii^ On 
the other hand, the mismatch problem is mitigated from the outset if contact layers made up of magnetic semicon- 

] ductors are usedi^i^iii Promisingly high values of the spin injection efficiency were obtained in experimentsi^iiSii^ 
• • ' using paramagnetic (II,Mn)VI DMS contacts. In the latter, the conduction band is spin-split due to the giant Zeeman 
effec1j2ii2Sii2^ in an external magnetic field, thereby making possible contact spin polarizations close to 100%. For this 
kind of spin injector, a novel magnetoresistance effect was observe d^^i^^ and theoretically analyzedi^i^liS^ within the 
^ ' standard drift-diffusion theory of electron transport. 

Spin-polarized transport in semiconducting structures formed of NMS layers and paramagnetic (II,Mn)VI DMS 
layers means electron transmission across spin-dependent potential profiles. These profiles depend on position via the 
position dependence of the electrostatic potential (composed of the conduction band edge potential and the external 
potential) and of the magnetic-field-induced spin splitting of the conduction band, which changes abruptly at the 
NMS/DMS interfaces. Transport across spin-dependent potential profiles corresponding to specific combinations of 
NMS and DMS layers has been studied theoretically in a number of case o^^i^^'^^ iS" in the (quantum-coherent) ballistic 
limit. Recently, the transport properties of II- VI resonant tunneling devices coupled to paramagnetic DMS contacts 
were found^i^ to depend strongly on the magnitude of the applied magnetic field. 

In the present work, we consider arbitrarily shaped potential profiles including internal and external electrostatic 
potentials and exhibiting arbitrary, position-dependent spin splitting, thereby covering any combination of NMS 
layers and paramagnetic DMS layers that may occur in semiconducting structures. We study spin-polarized electron 
transport within the (semiclassical) thermoballistic description'^'^ of carrier transport in nondegenerate semiconductors. 
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The basic element of this description is the thermobaUistic current, in which electrons move ballistically across intervals 
enclosed between points of local thermal equilibrium. The contribution to the current of each such "ballistic interval" 
is governed by the momentum relaxation length ("mean free path"). The thermobaUistic transport mechanism is 
intermediate between the diffusive (Drude) process, where the electrons move from one state of equilibrium to another, 
infinitesimally close- lying state (mean free path tending toward zero), and, on the other hand, the ballistic, collision- 
free motion across the electrostatic potential profile (mean free path tending toward infinity). In the thermobaUistic 
current, diffusive and ballistic transport are thus unified in a way that allows the effect of these two aspects of the 
transport mechanism to be studied. In Ref. [3^ . we have introduced spin relaxation into the thermobaUistic description 
for the case of electron transport in NMS. There, we have disregarded any spin splitting of the conduction band. 
Applications^i^^ have dealt with spin-polarized transport in heterostructures formed of an NMS layer sandwiched 
between two ferromagnetic metal contacts. In the present paper, we put forward the systematic extension of the 
thermobaUistic approach to spin-polarized electron transport across spin-split potential profiles. While we are mainly 
concerned here with the general formulation of this extension, we also treat, within a simplified picture, spin-polarized 
transport in heterostructures formed of an NMS layer enclosed between two DMS layers. Applications to specific 
experimental configurations will be deferred to future work. 

As a prerequisite to the thermobaUistic description, we formulate, in the next two sections, the details of ballistic 
spin-polarized electron transport across spin-split potential profiles. In Sec. II, we introduce the electron densities 
at the points of local thermal equilibrium enclosing a ballistic interval, and construct the currents injected from 
these points into that interval. In Sec. Ill, we continue the injected currents and the densities from the points of local 
thermal equilibrium into the ballistic interval, and introduce the balance equation that describes spin relaxation inside 
the interval. The off-equilibrium ballistic spin-polarized current and density are expressed in terms of a spin transport 
function that is related to the splitting of the spin-resolved chemical potentials. In Sec. IV, the thermobaUistic currents 
and densities are constructed by summing up the contributions from all ballistic intervals. The magnitude of these 
contributions is governed by the momentum relaxation length. The current and density spin polarizations as well as 
the magnetoresistance are each expressed as the sum of an equilibrium term determined by the spin-relaxed chemical 
potential, and an off-equilibrium contribution given in terms of the spin transport function. The spin-relaxed chemical 
potential is described through a resistance function. The latter, as well as the spin transport function, are each 
calculated from an integral equation. Section V deals with the thermobaUistic description of spin-polarized electron 
transport in DMS/NMS/DMS heterostructures. We derive explicit expressions for the current spin polarization and 
the magnetoresistance and compare these, in the limit of small momentum relaxation length, to those of the standard 
drift-diffusion approach to electron transport. Numerical results for DMS/NMS/DMS heterostructures are presented 
and discussed in Sec. VI. In Sec. VII, we summarize the contents of this paper and make some concluding remarks. 



II. POINTS OF THERMAL EQUILIBRIUM: DENSITIES AND INJECTED CURRENTS 

The thermobaUistic description of spin-polarized electron transport in NMS^i makes use, in a one-dimensional 
geometry, of ballistic electron currents and densities in "ballistic intervals" between two points of local thermal 
equilibrium, x' and x" (xi < x' < x" < X2)', see Fig. [TJ The ballistic electron transport across such an interval 
is determined by the densities at the points of local thermal equilibrium and by the shape of the potential profile 
inside the interval. For the purpose of generalizing the thermobaUistic approach to spin-polarized electron transport 
in paramagnetic DMS, we introduce the thermal-equilibrium densities nu(a;') and ri|_|(a;") for spin-up (|) and spin- 
down (I) conduction band states, as well as the spin-dependent potential profiles 

e^,ax)=E,{x)±A{x)/2 (1) 

at position x G [a;i,X2]- Here, the (spin-independent) potential Ec{x) comprises the conduction band edge potential 
and the external electrostatic potential, and A (a:) is the Zeeman splitting of the conduction band due to an external 
magnetic field^^ [we restrict ourselves to considering a single Landau level whose energy is assumed to be included in 
Ec{x)]. In developing our formaUsm, we assume both Ec{x) and A(a::), and hence e-i^i{x), to be continuous functions 
of X in the interval [xi,a;2]. Abrupt changes in one or the other of these functions, which occur at the interfaces in 
heterostructures, may be described, in a simplified picture, in terms of discontinuous functions. This aspect will be 
illustrated in Sec. V below for the case of DMS/NMS/DMS heterostructures. 

In the thermobaUistic transport mechanism, momentum relaxation takes place exclusively and instantaneously at 
the points of local thermal equilibrium. Spin relaxation via spin-flip scattering processes, on the other hand, occurs 
during the electron motion across the ballistic intervals (since the electrons spend only an infinitesimally short time 
span at the points of local thermal equilibrium, it is only inside the ballistic intervals that they can experience spin 
relaxation). This separation of momentum and spin relaxation is in accordance with the D'yakonov-Perel' relaxation 




FIG. 1: Schematic diagram showing, in a one-dimensional geometry, a paramagnetic semiconducting sample of length S enclosed 
between two nonmagnetic metal contacts. The spin splitting A(a;) of the potential profile Ec{x) gives rise to the spin-dependent 
potential profiles e|_|(a;). The coordinates x' and x" denote points of local thermal equilibrium. Electrons injected at x'{x") 
toward the right (left) move ballistically across the profiles ei^i{x) to reach the end-point x"{x') of the ballistic transport 
interval The thermoballistic equilibrium electron currents J{x) and J-{x) are constructed from the common spin- 

relaxed chemical potential ii{x') [with boundary values /ii,2 at the contact/semiconductor interfaces] according to Eq. (|77|) . 
Correspondingly, the off-equilibrium thermoballistic current J_ (x) is constructed from the spin transport function A{x') [with 
boundary values Ai^2]- 



mechanism^i^, which is known^i^ to be the dominant relaxation mechanism at large donor doping levels and at high 
temperatures. Therefore, our approach covers a broad range of cases interesting from the point of view of experiment. 



A. Electron densities 



The electron motion in a ballistic interval is activated at a point of local thermal equilibrium at one end and 
terminated at another such point at the other end. In the semiclassical model, the points of local thermal equilibrium 
are characterized by a local chemical potential (" quasi- Fermi level"). In a nondegenerate system, the spin-resolved 
electron densities, n|_|(a;'), are given by 

nu(-') = ^e-''[^-^(^')-''TaM], (2) 

where fj,i^^i{x') are the spin-resolved chemical potentials, Nc/2 ~ (27rm*//3/i^)'^/^ is the effective density of states of 
either spin at the conduction band edge, m* is the effective mass of the electrons (which, for simplicity, is assumed 
to be independent of position and of the external magnetic field), and /3 — {kBT)~^. In the spin-relaxed state, 
the potentials /i|^|(x') become equal to the common " spin- relaxed chemical potential" jj.{x'). Thus, the spin-relaxed 
densities at the point of thermal equilibrium x' , denoted by ri|^|(x'), are 
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For the Boltzmann factors e ^'^T,i(^ )^ we introduce the notation 

d^,^{x')^e-P^^'^^^'^ (4) 

and define 

d±{x')^d^{x')±d^{x') (5) 

(the "±" notation will be used mutatis mutandis for various other quantities appearing further below), so that we 
have 

n±{x') = ^d±{x')eP~^^^'\ (6) 

The total (i.e., spin-summed) spin-relaxed density, h{x') = n+(a;'), is 

h{x')^^d+{x')e^~^(-'\ (7) 

and for the "static" spin polarization P{x') of the conduction band electrons, i.e., the polarization in the spin-relaxed 
state at zero external electric field, we have 

We note the relation 

d^,,{x')^^[l±P{x')]^^{x') (9) 

expressing the Boltzmann factors 'di.i{x') in terms of P{x'). 

For the dynamical description of spin-polarized transport, we now introduce the ratios a-f4(a;') of the densities ^ 
to their spin-relaxed counterparts ([31), 

« (^') ^ = e'5[^T,i(-')-A(-')] . (10) 

working with these quantities will provide us with a description in terms of linear equations, instead of the nonlinear 
description that we would encounter when employing the chemical potentials /i|^|(a;') themselves. This aspect has 
been emphasized previouslyi^ within the diffusive approach to spin polarization in nondegenerate semiconductors. 
[We note that the quantity a-f j(x') introduced in Eq. pU)) is different from that of Eq. (3.8) in Ref. [s^; in the spin- 
relaxed state, we have in the present case au(a;') — 1, while in Ref. IH we have a^^i{x') — 1/2.] Using Eqs. (|4]), (O, 
(ini), and (Uni), we can write the densities n|.|(x') in the form 

n^„l{x') = i nix') [1 ± P{x')] aT,i(x') , (11) 

so that 

n±(x') = i h{x') [a±{x') + P{x') a^{x')] . (12) 

Here, the functions a^(x') and a-{x') appear as independent dynamical quantities. 

We now require the total density, n{x') = n^(x'), to be unaffected by spin relaxation, i.e., we set 

n{x') = il{x') . (13) 

This requirement is justifiedi^ for n-doped (unipolar) systems, for which spin-flip scattering processes involving valence 
band states can be disregarded. Substituting Eq. p3p in the upper Eq. (I12p . we find a+(x') and a^{x') related by 

a+{x')^2-P{x')a^{x'), (14) 

so that we can eliminate a+(a;') from the lower Eq. (jl2[) to obtain the spin-polarized density, n_(a;'), as 

n^ix') = nix') {Pix') + i[Q(a;')]'a-(x')} , (15) 
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where 

In the spin-relaxed state, we have 



"Ta(^') = i, (17) 

and hence 

a_(x') =0, a+{x') = 2. (18) 

Then, from Eq. p5p . the equiUbrium spin-polarized density at the point of thermal equilibrium x' , n^{x'), is 

n-{x') ^ n{x') P{x') (19) 

[see Eq. (j5])], and the off-equilibrium spin-polarized density, n_(a;'), is there 

h^{x') = n^{x') - n^{x') = i?i(2:') [Q(x')]^a-(a;') . (20) 

The off-equilibrium spin-polarized density is proportional to the quantity a_(x'), i.e., to the difference of the spin- 
resolved densities relative to their spin-relaxed values [see Eq. ([TU)) ]. but also contains, via the factor [(3(a;')]^, the 
effect of the static polarization P{x'). 

For later use, we establish the relation between the difference of the spin-resolved chemical potentials, /i_(x'), 
and the quantity a_(a;') which will appear as the key element determining the spin dynamics of off-equilibrium 
spin-polarized transport in the thermoballistic approach (see Sees. III.B and IV. C). From Eq. ([T^. we obtain, using 
Eqs. (0), (H), and fTO]), 

d^{x')a^{x') -h ^i{x')a^{x') = d+{x') , (21) 

and hence, using Eq. (|9|). 

a^^l{x') = l±\[lTP{x')]a,{x'). (22) 

With the help of Eq. (fTO]) , we then have 

1 (a^{x')\ 1 ( l + \[l-P{x')]a.{x') , 
'-^^ ^^-P'-" l^j = -P ll-i[l + P(.')]-(-') ' ' ^''^ 

and, reversely. 



« (^') = 2tanh(;3M-(^0/2) 

^ l + P(x')tanh(/3^_(a;')/2) ■ ^ ^ 

For the mean chemical potential jl{x') = ^fi-^-{x'), we find 
fl{x') = fi{x') + ^\n{aT{x')ai{x')) 

^ i,{x') + ^ ln({l + i[l - P(.t')] «_(a:')}{l - ^[1 + P(x')] a_(x')}) 

= /^(•^') + ^ 1^(1 - ^(2^')«-(^') - \[Q{x')n^-{x')r) . (25) 

Finally, we quote the relation between the spin-resolved chemical potentials fi^^i{x') and the common spin-relaxed 
chemical potential il{x'), 

i [1 + Pix')]e^''^ (-') + 1 [1 _ p(a;')]e^'^i = e'^'^^^'^ , (26) 
which, like Eq. ([21]), follows from Eq. but this time using Eqs. ([2]), dH), and JH). 
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B. Current injection 



In the thcrmoballistic description, a point of thermal equihbrium hes between two baUistic intervals (except at the 
ends of the sample, which will be considered further below), and electrons from the left- and right-lying intervals enter 
into it to be equilibrated instantaneously. At the same time, equilibrated electrons, which have no preferred direction 
of motion, are injected symmetrically (i.e., "half-and-half) into either ballistic interval, forming the spin-resolved 
electron current densities (electron currents, for short) Vgnijj^(x'), where — (27rm*/3)~^/^ is the emission velocity. 
If, however, there lies a potential barrier e^^{x' , x") — e|.|(a;') > (e™^{x' , x") is the maximum of the spin-dependent 
potential profile in the interval [x' , x"]) between, say, the left point of current injection, x' , and the opposite point of 
current absorption, x" , part of the current will be reflected at this barrier (and equilibrated while it travels back and 
forth between x' and the left side of the barrier, - we may call it the "confined current"). The part that surmounts 
the barrier (henceforth called the "injected current" ) is given by 



AT 

v,n^Jx')Tl^(x',x"), (27) 



J\jx',x") = ^(2/3/7rm*)i/2 / dp ^ e-'3[p'/2"*+^T,i(-')-MT.i(-')] e(p2/2„j* + „ (^'^ ^")) 

' 2 Jr. m 



where 



ri 1 (x', x") = e-'3[4':i(^'.^")-^Ta(=^')] ^ ) (28) 

i?t,x(a;') 

is the (classical) transmission probability for the spin-resolved current injected at the left end-point x' of the ballistic 
interval to reach the opposite right end-point x" . Here, we have introduced the notation 

i9'^^^{x',x") = e-f^'T.i(^'-^") , (29) 

which corresponds to the definition ^ for the Boltzmann factors i?|,|(a;'), with e|^|(a;') replaced with £™^{x',x") 
[note that the factors z9™(a;',x") and ■i9™(x',x") are local functions depending, for a given ballistic interval [a;',a;"], 
each on the position of the maximum of the potential profiles ei{x) and ei{x), respectively]. If the potential profile 
is constant along the interval or if its maximum lies at the injection point x' itself, then e^^(x' , x") = e|.|(a;') and 

T[^(x',x") = l. 

In analogy to Eq. pip for the spin-resolved densities, we now write the currents (|27p in the form 

J\jx',x") = lj'ix',x")[l±P"\x\x")]a^^^ix'), (30) 

where 

J'(x', x") EE j;(x', x") = ^^^{x', x")e'^'^(^') (31) 
is the total spin-relaxed current injected at x' into the interval [a;', a;"] from the left, and 

P-(x-,x")= ; . (32) 
wlj. (x', x") 

is a static, "nonlocal" spin polarization depending on the (generally different) positions of the maximum of the 
potential profile for the spin-up state and for the spin-down state, respectively [compare with the analogous relation 
([8]) for the local spin polarization P{x')]. 
From expression (|30p. we obtain, in analogy to Eq. (I12p . 

4(x',x") = ^j\x\x")[a±ix') + P'^{x',x")a^ix')]. (33) 

In line with the condition (|13p imposed on the total density n{x'), we now require the total current J^{x',x") = 
) to equal its spin-relaxed limit, 

f{x',x")^j\x',x"). (34) 
Following the argument leading to expression psp for the spin-polarized density n^{x'), we then find 

Ji{x',x") = f{x',x"){P"\x',x") + ^[Q"'{x',x")]^a-ix')}, (35) 
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where 

(see the analogous relation (fTH]) for [Q(a;')]^). 
In the spin-relaxed state [see Eq. (|18p ]. Eq. p3p yields for the injected equilibrium spin-polarized current 

JL{x',x") = f{x',x")P"^x',x"), (37) 

so that the injected ofT-equilibrium spin-polarized current is 

JL{x',x") = JL{x',x") - JL{x',x") = ij'(x',x") [Q"(x',x")]'a_(x'), (38) 

in parallel to Eqs. (fTO]) and ([20]) . respectively. 

All the preceding formulas hold also in the case where the left end-point of the ballistic interval coincides with the 
left end-point of the sample, x' = xi, and similarly for the right end-points, x" — X2, with the understanding that 
the chemical potentials il{xi) = jli and fi{x2) = jjL2 have fixed values given by the chemical potentials in the contacts 
bordering on the semiconducting sample at either end. 

III. SPIN-POLARIZED TRANSPORT ACROSS A BALLISTIC INTERVAL 

The electron currents discussed up to now are the currents injected into the ballistic interval at the point of thermal 
equilibrium as they leave the latter. Once they have entered this interval, they become transmitted ballistic currents, 
which will be considered in the following. 

A. Continuing the injected currents and densities into the ballistic interval 

We again consider injection at the left point of thermal equilibrium, x', of the ballistic interval [a;',a;"]. If the 
system is in the spin-relaxed state at the injection point a;', we have /i|_|(a;') = il{x') and a^^i{x') = 1. Then the 
transmitted (in the following, we omit the attribute "transmitted") ballistic spin- resolved current at x\ given by 
expression (|30p . propagates along the ballistic interval without spin relaxation, i.e., it is conserved. If, however, 
MT,i(^') 7^ i^{x'), the injected current will relax along the ballistic interval and hence will no longer be conserved. It 
is then natural to continue expression pop into the ballistic interval [x' , x"] by replacing a^_i{x') with a more general 
quantity ^(a;', a;"; x) depending both on the end-points x' and x" and on the position x £ [a;',x"], so that the 
ballistic spin-resolved current is obtained as 

J\^^{x',x"-x) = \ f(x\x")[l± P'^ix' ,x")]a\^^{x' ,x";x) = ^?^^(x', x") e'^'^^^'^ (a;', x"; x) , (39) 

with the initial condition 

a\^^{x',x"]x')^a^^l{x'), (40) 

which holds for all a;' > a:i. 

In line with the density ratios Q!| |(a;') in Eq. (|30p . the generalized quantities |(a;', a;"; a;) express the deviation of 
the currents inside the ballistic interval from their spin-relaxed values, so that here the condition for spin equilibrium 
is again 

a^-^ ^^{x , x" ] x) — 1 (41) 

[see Eq. (IT7|) ] , and hence 

a'_(a;',a:";a;) = , a\{x' ,x" ]x) ^ 2 , (42) 

as in Eq. ([T5)) . 

Extending condition (|34p and the procedure ensuing therefrom, we obtain for the (conserved) total ballistic current 
inside the ballistic interval 

J'(x',a:";a:) = f_^{x',x"-x) = j\x',x") = !^^!^(x', a:")e^'^("') (43) 



8 



and, similarly, for the (conserved) ballistic equilibrium spin-polarized current 



f_{x',x";x) = j'(x',x")P"(x',x") = !^^-^{x',x")e''''^^'^ (44) 



[see Eq. (I37p ]. while for the ballistic off-equilibrium spin-polarized current we have 

J'_ix',x";x) = y\x',x")[Q"'{x\x")fai{x',x";x) 

^i9"(x', x")ef^i'^^'^ai{x', x";x) (45) 



2' 
VeNc 



[see Eq. dSHl)] , where 

2i9™(x',x")i9r(x',2;") 

'^'"("''^"^ ^ \^{x',xn ^ hKi-\^")[Q"'{:>:\^")r- (46) 



In Eq. (|45p . the spin dynamics of the ballistic off-equilibrium spin-polarized current j'_{x' ,x";x) is determined by 
the factor aL{x' , x"; x), whose dependence on the coordinate x reflects the spin relaxation inside the ballistic interval 
[x',x"]. This dependence will be considered in Sec. III.B. 

The density ^{x',x";x) associated with the ballistic spin-resolved current Jj ^{x',x";x) will be called "ballistic 
spin-resolved density". It is obtained from Eq. (l39|) by replacing Ve with C'^^i{x' ,x";x)/2, 



J 



nl^{x',x";x) = -^D'(^^{x',x";x)e'^^'^'''^ a\jx',x";x) (47) 

(see the analogous relation between the current (2.2) and the density (2.7) in Ref . [33) . Here, 

D^l^ix',x";x) = C^,^{x',x";x)r^]^ix',x") (48) 

and 

CT,i(^',^";^) . eM[-\nT,A^',x";x)]'/^) (49) 

T],i{x',x";x) 



[see Eq. (2.8) of Ref.[3|, with 



,9m ( f ii\ 



(50) 



the latter quantity may be interpreted as the transmission probability corresponding to injection at the point a; toward 
the region containing the maximmn of the potential profile in the interval [x' , x"]. The ballistic velocity v^^^l{x\ x" \ x) 
is given by 

v^,l{x',x"-x) ^ \ \ ^ „ 51 

nY^{x',x"\x) C-\,i{x' ,x";x) 

It is not affected by spin relaxation, since the spin-flip mechanism is assumed not to influence the kinematics of the 
electron motion. 

From Eq. (|47|) . the continuation of half (see the remarks at the beginning of Sec. II. B) the densities (fT2|) into the 
ballistic interval now follows as 

n'±(a;', x"; x) = \n\x' , x" \ a;)[a'±(a;', x"; x) + P^{x\ x" \ x) a^(a;', x" \ x)\ , (52) 

where 

-III II •■ -1,1 If ^ '^r. 



n 



{x',x";x) = n^+{x',x";x) = -^D"\x' ,x";x) 6^^^"='^ (53) 



[see Eq. ([7])] is the total ballistic spin-relaxed density, and, in generalization of expression for P"''{x\x") 
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Extending condition (jl3p into the ballistic interval by requiring 

1 1 I II \ I I I II \ -1/ I II \ /'c;t;^ 

n[x,x ■,x) = n^[x,x ■,x) — n[x,x ;x), (55) 

we find from Eqs. ([5^ 

a'^{x',x";x) ^2- P^{x',x";x)a'_{x\x";x) (56) 
[see Eq. ([T4|) ]. and hence for the ballistic spin-polarized density inside the interval [x' , x"], in parallel to Eq. p^ . 

nL(x', x";x) = n'ix\ x" ; x){P'^\x' , x" ; x) + \[Qc{x\ x" , x)fa^^{x\ x" ■ x)} , (57) 

where 

, ADfix',x"\x)DV''{x',x":x) 
[Q^{x\ x", x)f ^ 1 - {P-S{^'. - ^ ' ,/ \i (58) 

[yj™ (x', a;"; x)Y 

[see Eq. p6p ]. Thus, we obtain for the ballistic equilibrium spin-polarized density 



hi {x\ x": x) = -^D"^{x', x"; x)e^^^^"> (59) 



[see Eq. (|19p ]. and for the ballistic off-equilibrium spin-polarized density 

N. 



Ux\x";x) = -^D"\x\x";x) e^'^^^'^aUx' ,x"-x) (60) 



[see Eq. ^^], where 



, „ 2DV^(x',x";x)DY'(x',x";x) , , „ , „ o 

D^{x',x";x) = ' ] y ' = \DX{x\x"-x)[Q^{xix"-x)]\ (61) 

D_l[x' ,x";x) 



in parallel to Eq. ((46l) . 

We remark that for constant potential profiles, e|_|(a;) = e^ i, when C^^iix' ,x";x) = 1 so that D^{x',x";x) = 
•d™{x' ,x"), Eq. ([55]) becomes •n}{x' ,x" \x) = h{x')/2. This reflects the fact that the left-hand side of this relation 
refers to the density associated with the "half-sided" injected current, while h{x') is the total density at the point of 
thermal equilibrium x' . 



B. Spatial behavior of the off-equilibrium spin-polarized current and density 

The total ballistic current J'(x',a;") and density ti''{x', well as the ballistic equilibrium spin-polarized 

current Ji(a;', x"; x) and density ni{x' , x"; x), are determined by the quantities 'd''J!:{x' , x") and D'!^{x' , x"; x) [which, 
in turn, are completely determined by the potential profiles e-f^|(x')], and, most importantly, by the spin-relaxed 
chemical potential /i(x'), which is the only dynamical quantity appearing. The latter is to be calculated according to 
the thermoballistic procedure presented in Ref . "s^ (see below) . 

By contrast, the ballistic off-equilibrium spin-polarized current ,r_{x' ,x" \x) and density n'_(x', x"; x) contain, in 
addition, the common factor aL(x',x";x), whose spatial behavior (x-dependence) is determined by the process of 
spin relaxation in the ballistic interval [x',x"]. This process is governed by the balance equation connecting the 
off-equilibrium ballistic spin-polarized current and density (see Ref. 34), 

^Jl{x',x"-x)+'^hi{x',x"-x)=Q- (62) 
dx Is 

here, Is = 2veTs is the (ballistic) spin relaxation length, which we assume, as in the case of the effective electron 
mass m*, to be independent of position and of the external magnetic field. The spin relaxation time is given by 
1/ts = l/'T'Ti + l/'T'iTj where l/r-fj^ (l/'''!?) the rate for spin-flip scattering from spin-up (spin-down) to spin-down 
(spin-up) states. Inserting expressions ([151) (|60p in Eq. ((5^ . we now obtain a first-order differential equation for 
aL(x', x"; x), 

-^aiix',x";x) + ^^^^^^fi^a'_ix',x";x)^0, (63) 
ax Is 
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where 

rir I II \ — ri f I ii \n r i n \ "^^{^ ) [Q™ (^'' ■^)]^ D^{x,X ; x) 

Clx ,x ■x) = Cf{x,x ^xjCiix ,x ;x)——f— — - — - — — — — — — J — ■ (64) 

^ ' ' ^ IV : , / u > > ' D'^ix^x'^ix) [Q"'{x',x")]^ ^ ' 

The solution of Eq. ((63|) obeying the initial condition pO)) is 

[x', x"; x) = a_ {x') e-c(-'.-";-'.-)/'s , (65) 

where 



C{x,x";zi,Z2)= dzCix',x";z), (66) 



with z< = min(zi,Z2), 2:> = max(zi,Z2). 

For the ballistic off-equilibrium spin-polarized current, we then have from Eq. (|45 



Ji(x',x";x) = !^^"(a;',a;")A(a:')e-'^("''""'"''"'/'% (67) 



and for the ballistic off-equilibrium spin-polarized density from Eq. ((60 



hlix\x";x) = ^i?'"(a;',x";x)^(a;')e"'^^"''""'"''"^/''- (68) 

Here, we have introduced the spin transport function 

A{x') = e'^'^(="') a_(x') ^ e''''T(^') _ ef^mi^') _ (69) 

[Owing to the different definitions of a^,i{x') in the present work and in Rcf. 34, the quantities A{x') are not the same 
in the two papers.] In the thermoballistic approach, this function completely describes the dynamics of off-equilibrium 
spin-polarized transport. In terms of the difference of the spin-resolved chemical potentials, fi-{x'), which commonly 
has been used as the basic dynamic variable in previous approaches to spin-polarized transport, it is expressed, owing 
to Eq. dUl), as 

^ 2eM-'hanh{/3^i^{x')/2) 
^ ' l + P(a;')tanh(/3^_(x')/2) ■ ^ ' 

The procedure for the calculation of the spin transport function will be outlined in Sec. IV. C. 

The spatial behavior of the off-equilibrium spin-polarized current and density in the ballistic interval i.e., 
their dependence on the coordinate x, is determined by the factor e~^'^^ ■.x)lis ^ Hence, spin relaxation in this 
interval departs from a purely exponential behavior unless the potential profiles t\^i(x) are constant over the interval. 



C. Net ballistic currents and joint ballistic densities 

So far, we have only considered injection from the left end-point, a;', of the ballistic interval [x', x"\. The expressions 
for the currents and densities obtained for that case are easily transcribed to the analogous expressions for injection 
from the right end-point a;", J''(x',x"), J^(x' ,x";x), J!l(a;', x"; x), n^{x',x";x), rV^{x' ,x";x), and ff_{x' , x" ] x) , by 
replacing everywhere the superscript / with r, and interchanging the arguments x' and x" [owing to the symmetry of 
e|"|(x', x"), the quantities d^{x\ x"), '&"^{x' , x"), D™{x' , x"; x), and D™{x' , x"; x) are symmetric under the exchange 
of x' and x"] . 

Now, combining the ballistic currents injected from the left and right, we have for the (conserved) net ioia/ ballistic 
current inside the ballistic interval, using Eq. (|^ . 

J{x',x") = J\x',x") - r{x',x") ^^-^d^{x',x") [e'''^("') - e'5'^("")] . (71) 
Using Eq. ()44|) , wc find for the net ballistic equilibrium spin- polarized current 

J-{x\x") = JL{x' ,x") — JL{x' ,x") = J{x' ,x") P^"'{x' ,x") 
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while for the net baUistic off- equilibrium spin-polarized current, we have, using Eq. (j67p . 

J — (^X ^ X ^ X^ - — J _ ^ X ^ X J 3j ^ " J (yXj ^ X ^ X ^ 

Similarly, using Eq. ([55]) . we find for the joint total ballistic density associated with the net total ballistic current ([7T|) 

n(x', x"; x) = n\x', x"; x) + n'-(a;', x"] x) = ^D"^{x', x"; x) [e'^'^^"') + e>^^'^^"^] , (74) 

while for the joint ballistic equilibrium spin-polarized density we obtain, using Eq. (|59p . 

n_(x', x"; x) = rlL{x', x"; x) + nL [x' , x" ; x) = n{x\ x" \ x) P^{x' , x" \ x) 

= ^i:)™(a;'^a;";a;)[e'3A'(^')+e''/i(:^")]. (75) 

Finally, for the joint ballistic off -equilibrium spin-polarized density, we have, using Eq. (|68p . 
n-{x',x"]x) — n''_{x' , x" ; x) + fL^_{x' , x" ; x) 

= ^Z?™(x',x'';a;)[A(a;')e-'^(^'^^'''^''^)/'= +^(a;'')e^'^(^''="'''^'^'')/'=]. (76) 

It becomes apparent from the preceding results that the role of the x-independent, local Boltzmann factors i?™|(x', x") 
in the currents is taken over, in the densities, by the x-dependent, nonlocal factors DJ^^{x',x";x). 

IV. THERMOBALLISTIC SPIN-POLARIZED TRANSPORT 

The thermoballistic approach to semiclassical electron transport has been developed in Ref. "sS, and generalized in 
Ref. [34: to describe spin-polarized electron transport across a spin-degenerate conduction band. In this approach, the 
net ballistic currents and joint ballistic densities defined in ballistic transport intervals form the building blocks for 
the construction of the corresponding thermoballistic currents and densities. In order to deal with the case of a spin- 
split conduction band, we have to modify and augment the formulation of Ref. Isil. We construct the thermoballistic 
currents and densities in terms of the net ballistic currents and joint ballistic densities derived in Sec. III.C, and 
outline the procedures for calculating the equilibrium chemical potential and the spin transport function. From 
the thermoballistic currents and densities, the corresponding equilibrium and off-equilibrium spin polarizations are 
determined, and the magnetoresistance is obtained. 

A. Thermoballistic currents and densities 

Following Refs. Issl and [3^. we introduce, for given position x inside the semiconducting sample extending from xi to 
X2 (see Fig.[T]), the set of ballistic intervals [x' , x"] that enclose x, thus restricting x' and x" by xi < x' < x < x" < X2- 
From the net ballistic currents and joint ballistic densities of Eqs. ([7T|) - ([75)) . here summarily represented by the 
symbol J-(x',x";x), the corresponding thermoballistic currents and densities, J-{x), are constructed by summing up 
the currents and densities !F{x\ x"; x) over all ballistic intervals. In the summation, the ballistic intervals are weighted 
(assuming one-dimensional transport) with the probability e~'^ ~^ that the electrons traverse the interval without 
impurity or phonon scattering, multiplied by the probability dx' /I {dx" /I) of their being absorbed (or emitted) in an 
interval dx' {dx") at an end-point x' {x"). At the ends of the semiconducting sample, xi^2, absorption and emission 
occur with unit probability. Here, the quantity I is the momentum relaxation length (mean free path), which governs 
the relative contribution of the ballistic and diffusive ("thermal") transport mechanisms (see Ref. [33l ). Like the 
effective electron mass m* and the ballistic spin relaxation length Ig, the momentum relaxation length I is assumed 
to be independent of position and of the external magnetic field. 

By this procedure, the thermoballistic currents and densities J^{x) are obtained in the form 

T{x) = e-("^-"i)/'.F(a:i,a;2;a;)+ / — e-^^^"" .F(x', X2; x) 

J xi ' 

+ / '^_e-i-"-^)nT{x,,x";x)^^ ^e-^^-^^/'j^{x',x";x). (77) 

J X J xi ^ J X ^ 
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Thus, the total thermobaUistic current T{x) = J{x) is given by Eq. (|77|) with T{x' ,x";x) replaced with J{x',x") of 
Eq. ((7T|) . Comparin g w ith Eq. (2.12) of Ref. [S^, we see that the current J{x) of the present work corresponds to the 
current J{x) of Ref. [3J if there one sets = and substitutes for the weights Wn{x' , x"; I) [n = 1, 2, 3] the function 

w+{x',x";l) = ie-l"'-""l/'^!;^(x',a;"), (78) 

which in the limit of zero spin splitting, ??"'(a;', x") = 2, reduces to the (n-independent) weight Wn{x',x";l) for 
one-dimensional transport. 

The equilibrium thermobaUistic spin-polarized current J^{x) = J-{x) follows from Eq. (|77p by identifying 
J-(x' , x"\ x) with J-{x',x") of Eq. ([72| . For the current J^{x), there is no explicit counterpart in the formulation of 
Ref. |3J, in which the equilibrium part of the thermobaUistic spin-polarized current merely has been parameterized in 
terms of an arbitrary constant a_ [see Eq. (3.27) of Ref. [Sj]. 

The off-equilibrium thermobaUistic spin-polarized current J^{x) = J_ (x) is obtained from Eq. (|77p by replacing 
J-{x' ,x";x) with J_{x' ,x";x) of Eq. (|73p . The resulting current J^{x) is seen to correspond to the current given by 
Eq. (3.28) of Ref. S if the function 

w{x', x"; /) = i e-l"'^""l/' d"\x', x") (79) 

is substituted for the weights Wn{x' , x" ; I) . Here again, in the limit of zero spin splitting, d™{x' ,x") = 1, the function 
w{x' ,x";l) agrees with the weight Wn(x' ^x" \l) for one-dimensional transport. 

We note that the total and equilibrium spin-polarized currents, J{x) and J^{x), both are linear functionals of the 
factors d'^{x' , x"), and hence of the Boltzmann factors corresponding to the maxima of the potential profiles e'\,i{x) 
in the ballistic interval [a;',a:"]. The off-equilibrium current J_(x), on the other hand, is, via its dependence on the 
factor 'd™{x' , x"), a nonlinear functional of the Boltzmann factors. 

The thermobaUistic densities n(a;), n-{x), and fi-{x) corresponding to the thermobaUistic currents J(x), J-{x), and 
J-{x), respectively, are evidently obtained by substituting the joint ballistic densities of Eqs. ([7i)) " ((75|) for !F{x' , x"; x) 
in Eq. (|77|) . Com par ing the thermobaUistic densities of the present work to the densities given by Eqs. (2.17), (3.29), 
and (3.30) of Ref. 34, we find that the remarks made above in conjunction with the thermobaUistic currents also apply 
here if for the weights given by Eq. (2.55) of Ref. 34 the functions w^{x' , x"; /) and w{x\ x" \ I) of Eq. ([75)) and (|79p . 
respectively, are used, with i?"'(a;', a;") and i?"'(a;', a;") replaced with D"^ (x' , x" ; x) and {x' , x" ; x) , respectively. 



B. Spin- relaxed chemical potential and equilibrium spin polarizations 



The net total ballistic and ballistic equilibrium spin-polarized currents, Eqs. (|7ip and (|72p . respectively, as well 
as the joint total ballistic and ballistic equilibrium spin-polarized densities, Eqs. ([7^ and ([75)) . respectively, and 
hence also the corresponding thermobaUistic currents and densities, are dynamically determined by the spin-relaxed 
chemical potential jl{x). The procedure for calculating the latter quantity is set up^ by employing the conditions 

1 

/ dxJ{x) = J, (80) 

X2 - xi Ja,^ 

where J is the total physical current, and 

J{x+) = Jix^) (81) 

for the total thermobaUistic current J{x) given by Eq. ([77)1 with J-'{x) = J{x). 

We now rewrite J{x), with the help of Eq. ([7T|) . in terms of the (spin-relaxed) "resistance function"— 

^(a;) = ^ [e^A(-i) _ e/3A(-)] . (82) 

Inserting the result in Eq. (jSOp and replacing the parameter X2 with the variable a;, we obtain a linear, inhomogeneous, 
Volterra-type integral equation for x(2;). 



X — Xi 



r dx' 

f+{x;xi;l)x{x) + — K+{x,x';xi;l)x{x') ^ , (83) 



where 



r dx' 

f+{x;xi;l) = u+{xi,x;l) + —u+{x',x;l) (84) 

Jxi ' 
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and 



with 



dx" 

K+{x,x';xi;l) — u+{x' ,x;l) — u+{xi,x';l) + j —— u+{x' ,x";l) , (85) 



u+{x',x";l) = j—w+{x',x";l), (86) 

where w-f-{x' , x"\ I) is given by Eq. (|78p . At x = xi, the resistance function x{^) is discontinuous, 

X(xi) = 0, (87) 

x(a;f) 



(note that we assume the potential profiles 614(2;) to be continuous in the interval [xi, X2\)- The function f+{x; xi; I) 
and the kernel K+{x,x';xi;l) are seen to depend on the profiles eij^i{x) solely via the quantity ■d^J^{x' ,x"). The 
solution of Eq. (|83p is obtained, in general, by numerical propagation. 

In order to construct a unique chemical potential fi{x), following Ref. [H, we have to solve Eq. ([55]) twice: first, with 
the original profiles e-f j(a;) to obtain the solution xi{x)i ^^^id second, with the spatially reversed profiles ^(x) — 
e|^l(xi + X2 — x) to obtain the solution x^(x). Defining 

X2ix) = xU^i + X2 - x) , (89) 

we introduce the function 

X-{x) = ai[xi(a;) - 5X1(2:2)] - a2[x2{x) - ^X2{xi)] (90) 

and the constant 

X = aiXi{x2) + a2X2{xi) . (91) 

The coefficients oi 2 are given2i by 

ai,2 = ^ , (92) 

with 

dx' 

ai= —{w+{xi,x';l)[x2{x')-X2{xi)]+w+{x',X2;l)x2ix')}, (93) 
/■^^ dx' 

02=/ -j- {w+ixi,x';l)xiix') +w+{x',X2;l)[xi{x') - Xi{x2)]} , (94) 

and a = ai+a2, SO that 01+0,2 = 1. For symmetric potential profiles, e| |(x) = e|.|(x), we have X2(2;) = Xi{xi+X2—x), 
and hence oi = 0,2 = 

In terms of X-ix) and Xi the thermoballistic spin-relaxed chemical potential jl{x) is given^^ by 

e^'^(^) =r/+-2^^^77_ (95) 
X 

(xi < X < X2), where 

= i(e'3'ii ±e'3'^^) (96) 

and /ii_2 = /i(a;i.2). The chemical potential /i(x) thus obtained in the thermoballistic approach characterizes the 
thermal-equilibrium attribute of the thermoballistic system, which co-exists with its ballistic attribute. It immediately 
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leads, via Eqs. ([7]) and ((T9l) . to the total spin-relaxed density n{x) and the equilibrium spin-polarized density n^{x), 
respectively. 

By analogy with Eq. (2.44) of Ref. [s^, we introduce the (dimensionless) reduced resistance 

X^^d^{x,,X2)x, (97) 
which fixes, for given total current J, via the current-voltage characteristic 

J = ^VeN, ^^ixi,X2) e^^"' 4 (1 - e-^-^) (98) 

X 

[see Eq. (2.42) of Ref. [11] the "spin-relaxed" voltage bias 

V^^l^ (99) 

e 

between the contacts bordering the semiconducting sample. In the zero-bias limit, x determines, via its dependence 
on the Zeeman splitting A(x), the magnetic-field-dependent equilibrium resistance 



eJ 



^ (100) 



where nf^ = iv^ e-''['^<=(^i)~'^il is the total (spin- relaxed) electron density at Xi for zero external magnetic field, i.e., 
for A(a:) = 0. [Note that, owing to its definition in terms of a current density, R has the dimension of an interface 

resistance.] It is notationally convenient, and physically meaningful, to introduce, instead of the density n^i \ the 
Sharvin interface conductance^^ 



g['^ = ^e\,n^^\ (101) 



so that 



Adding to the spin-relaxed bias V the off-equilibrium bias, one obtains the magnetoresistance proper (see Sec. IV. D). 

The total thermoballistic current J{x) is obtained in terms of /i(x) or, equivalently, X-(^) cind x by substituting 
expression ([7T|) in Eq. (|77p with !F{x) = J{x). In an analogous manner, the total thermoballistic density n{x) 
associated with the total thermoballistic current is obtained by referring to Eqs. (|7I)) and ([77)) . 

We define the equilibrium current spin polarization Pj{x) and the equilibrium density spin polarization Pn{x) in 
terms of the thermoballistic currents and densities, respectively, as 

Pj{^) = ^ (103) 

and 

P„(x) = ^ . (104) 
n(x) 

The thermoballistic currents and densities are, in general, not equal to the corresponding physical currents and 
densities. However, the relative spin content is the same in both, so that the "thermoballistic" ratio in Eqs. (|103p and 
(I104p . respectively, is equal to the corresponding "physical" ratio. 

In view of Eq. ([55]) and the first of Eqs. ([55]) . the term in brackets on the right-hand side of Eq. ([7T|) for the net total 
ballistic current J{x',x") is proportional to the total physical current J, since the function X-i^) is independent of 
J. From expression (|77p . the total thermoballistic current J{x) is then again found to be proportional to J, and the 
same holds for the equilibrium thermoballistic spin-polarized current J-{x) [see Eqs. (j72p and (|77p ]. The equilibrium 
current spin polarization Pj (x) is therefore independent of J (it is an intrinsic property determined exclusively by the 
spin splitting of the conduction band inside the semiconductor). The equilibrium density spin polarization Pn{x), on 
the other hand, is found to depend, in general, on the bias V, i.e., on J [see Eq. (2.53) of Ref. ^3^]. 
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C. Spin transport function and off-equilibrium spin polarizations 



The off-equilibrium thcrmoballistic spin-polarized current J_(a;) is obtained by substituting expression (|73p in 
Eq. ([77)1 with T{x) = J-{x), and analogously for the corresponding density (fZS)) . Following Ref. 34, we invoke the 
spin balance equation connecting these two quantities and arrive at an integral equation for the spin transport function 
A{x) of the form [see Eqs. (3.35)-(3.37) of Ref. [H 

dx' - 

W[xi,x-lls) Ai+W{x,X2]l,ls) A2-W{x-xuX2]l) A{x) + / —W{x',x-l,ls)A{x')^Q, (105) 

where 

Wix', x"; I, Is) = wix', x"; I) e-c(-'.-" (106) 

and 

dx' 

W{x;xi,X2;l) = w{xi,x;l) + w{x,x2;l) + / -—w{x',x;l). (107) 

Jxi ' 

Here, Ai,2 — A{xi^2) are the values of the spin transport function at the contact side of the interfaces with the 
left and right contacts, respectively; they are determined by the chemical potentials (2:1^2) at these positions 
[see Eq. dMl]. The exponent C {x' , x" ; x' , x") is given by Eq. and w{x',x";l) by Eq. Equation (fTUS)) is 

a linear, inhomogeneous, Fredholm-type integral equation for the spin transport function A{x). The corresponding 
homogeneous equation is solved by A{x) = only, so that the solution A{x) of Eq. (|105p for a;i < a; < X2 is linear 
and homogeneous in Ai and A2. The function A{x) is not, in general, continuous at the interfaces with the contacts, 
A{xf) Ai,A{x2) 7^ A{x2), as a consequence of analogous discontinuities of the spin-resolved chemical potentials 
f^hlix) ["Sharvin effect"]. In general, Eq. ()105p can be solved numerically by using matrix methods after discretization. 
The off-equilibrium current spin polarization is given by 

Pj{x)^'-^, (108) 
J(x) 

and similarly, the off-equilibrium density spin polarization, by 

Pnix)^^. (109) 

Here again, as in the case of the spin-relaxed polarizations, we calculate the polarizations in terms of thermoballistic 
currents and densities, respectively. 

These off-equilibrium polarizations are determined, via Eq. (|77p . by the current J^{x' ,x";x) and density 
'h_{x' ,x";x), Eqs. (j73p and (|76p . respectively, which are linearly connected with the spin transport function A(x), 
and therefore, linear in the boundary values A12, 

= Fi{x) Ai + F2ix) A2 , (110) 

with "formfactors" Fi{x) and F2{x). 

It remains to determine the quantities Ai^2, or rather ai^2 = ct-{xi,2) [see Eq. ([69)) ]. We do this by making use of 
the continuity of the total current spin polarization 

Pj{x) ^ Pjix) + Pj{x) , (111) 

closely following Sec. IV of Ref. At the interfaces with the semiconductor at 3: = a:i and x = X2, respectively, the 
current polarizations Pj{xi^2) in the left and right (semi-infinite) external contacts, which we take as nonmagnetic 
metals, are found as 



here. 



P.(x.2) = T^lnri±|ii^^:H]^); (112) 



L 
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where ai^2 are the conductivities of the two contacts, Ls ' ' the spin diffusion lengths, and Pi^2 = P{xi.2)- We 
note that Eq. (|112p differs from Eqs. (4.7) and (4.8) of Ref. [13 with respect to the role of the polarizations Pi^2' in 
the former equation. Pi 2 refer to the paramagnetic semiconductor, while in the latter equations they refer to the 
ferromagnetic leads (note also that the normalization of the quantities ai_2 used here differs from that used in Ref. 
The corresponding values of the polarization in the semiconductor are given by 

Pj (xia) = Pj {X1.2) + ^^[Fi{xi,2) Ai + F2{xi,2) A2] , (114) 

with Ki^2 — J{xi^2) I J ■ Equating the right-hand sides of Eqs. (jll2p and (jll4p . we obtain a pair of nonlinear equations 
which determine the quantities ai^2 = Ai^2e~^'^^^'-'^\ and thereby, the off-equilibrium current spin polarization Pj{x), 
in terms of the equilibrium polarization Pj{x), i.e., of the Zeeman splitting A(a;), and the material parameters of the 
metal contacts and the semiconductor. 

In the zero-bias limit, we have |ai.2| ^ 1, and the pair of equations for ai^2 become linear, 

T S4 «i,2 - ^^2^^^ [Fi(xi,2)ai + F2{xi.2)a2] = Pj(xi,2) J • (115) 

As solutions of Eq. (I115|) . the quantities ai.2 are proportional to the current J, so that the current spin polarization 
Pj{x) is independent of J. 

Summing up, we obtain the current spin polarization along the entire heterostructure, Pj{x), as follows. In the 
nonmagnetic contacts, it decays exponentially away from the interfaces at xi.2 with decay length Ls , 

^./(-) = Tg5A^-(x,2)e-l— 1/-^^'" (116) 

[see Eqs. (4.3) and (4.5) of Ref. 34 with Pi. 2 = 0]. In the semiconductor, it is given by expression piip . with the 
formfactors ^1.2(2^) determined by the solution of the integral equation (jlOSp . The explicit form of the latter depends 
on the structure of the semiconductor. 

D. Magnetoresistance 

Once the quantities ai^2 are determined, we can obtain the magnetoresistance R = [V/eJ]j^o, where V is the 
voltage bias related to the mean chemical potentials jli^2 = ^(3:1,2) at the contact side of the contact-semiconductor 
interfaces, 

V^-{fli-fi2) (117) 

e 

(here, we assume the metal contacts to have infinitely high conductivity, so that there is no voltage drop inside the 
contacts). From Eq. (|^ . we have 

Mi,2 = p-ia + ^ ln(a|(a;i,2)aj.(xi,2)) , (118) 

which in the zero-bias case reduces to 

A*l,2 = Mi,2 - 7^ A,2ai,2 ■ (119) 
zp 

Subtracting the two equations (|119p from one another, we obtain the magnetoresistance as 

R = R + R, (120) 
where the equilibrium contribution R is given by Eq. (jl02p . and 

R^-^^{Piai~P2a2) (121) 

is the off-equilibrium contribution. 
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V. DMS/NMS/DMS HETEROSTRUCTURES: THEORY 

As an illustrative example, we now consider the thermoballistic description of spin-polarized electron transport in a 
heterostructure formed of an NMS layer sandwiched between two DMS layers. As a structure of this kind is composed 
entirely of semiconducting material, one should ideally treat it as a single sample, so that the ballistic intervals [x' , x"] 
may enclose one or both of the DMS/NMS interfaces. In this case, one has to consider carefully the effect of the 
interfaces on the electron motion. The interfaces are characterized by (i) abrupt changes in the material parameters, 
reflecting changes in impurity and phonon scattering and in magnetic scattering, (ii) high structural disorder in their 
vicinity, and (iii) an abrupt change in the electrostatic potential due to the change in spin splitting. Taking into 
account the combined effect of these features in a consistent way is beyond the scope of the present work. 

Since, here, we place emphasis on the general formulation of the theory, we adopt a pragmatic point of view: we 
assume the different layers in a DMS/NMS/DMS heterostructure to be homogeneous and require the interfaces to act 
as fixed points of thermal equilibrium. That is, we apply the thermoballistic description separately to the different 
layers (thus dealing with conduction band potentials and spin splittings that are, in general, discontinuous at the 
interfaces) , evaluating for each layer the spin transport function and the current spin polarization for a homogeneous 
semiconductor, and subsequently connect these functions across the interfaces to obtain the position dependence of 
the spin polarization as well as the magnetoresistance. In a previous publication^^, we have applied this concept 
within a heuristic approach. 



A. Homogeneous semiconductor 

Inside a homogeneous semiconducting layer, at zero bias, we have a flat spin- independent potential profile, Ec{x) = 
Ec, and a constant spin splitting, A (a;) = A, so that the profiles e-\^i{x) are given by 

e^^^{x) = E^_±A/2 (122) 

{xi <x< X2)- For the functions d'^{x',x") and ?9™(x',a;"), respectively, we then have from Eq. ([5]) 

d"^{x\ x") = 2e~>^^^ cosh(/3A/2) = ^e'^'^^ (123) 

and 

op 

?9™(x', x") = -2e-^^= sinh(/3A/2) = -jre'^^^ , (124) 

where Q = (1 - p2)i/2^ ^nd P = -tanh(/3A/2) is the static spin polarization. Now, using Eqs. (|78|) and (I123P 
in Eq. we easily see from Eqs. that the resistance functions Xii^) ^^id X2{^) for the homogeneous 

semiconductor are given by the corresponding functions for Ec{x) = A(x) = [see Eq. (2.68) of Ref. IsJ], multiplied 
by e'^^^Q. Hence, the function X-{^)/Xh is independent oi Ec and A. Equations ([95|l and (|96|) show that, for given 
//i^2, the spin-relaxed chemical potential jlh{x) is also independent of Ec and A. Using Eqs. (|7ip and p23|) in Eq. (|77p . 
we then find the total thermoballistic current in the semiconductor as 

Jh{x) EE ,h - Jf ^ , (125) 

where jj^'' is the (conserved) total physical current for Ec ~ A — 0. Similarly, using Eqs. ((72|) and (|124p in Eq. (|77|) . 
we have for the thermoballistic equilibrium spin-polarized current 



j!!(x)^ J'! -e-'3^=^4"\ (126) 



and hence for the zero-bias equilibrium current spin polarization of the homogeneous semiconductor 



thereby retrieving the static spin polarization, as was to be expected. 
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As to the ofF-equilibrium current spin polarization, we have from Eq. (|46l 



this factor drops out from the integral equation (jlOSp . and as its solution, the spin transport function A{x) is seen to 
be equal to the function for i?c = A = 0, given by Eq. (3.56) of Ref. 34, 

A{x) = Cie-(^-^i)/^ + C2e-(^^-^)/^, (129) 
where the generalized spin diffusion length L is given by 

L = yiF, = , (130) 



111 , , 

T = T + L' (1^1) 



with 



and 

is = (132) 

is the spin diffusion length proper. The coefficients Ci^2 are linearly connected with the boundary values Ai^2 [see 
Eqs. (3.57) and (3.58) of Ref. '3^. 

Using Eq. (|73p with Eq. (|128p in Eq. ([TT]) . we find for the thermoballistic off-equilibrium spin-polarized current 

J'lix) ^ Qe-"^^ j'l^°\x) , (133) 

where J^*-*'-' (x) is given by 

rj'\x)^-v,NJ^ (134) 
ax 

[see Eq. (4.11) of Ref. [13; the factor 2 appearing in the right-hand side of the latter equation reflects the fact that 
the normalization of A{x) used there differs from that used in the present paper]. Combining expressions (|133p and 
(I125p . we obtain the zero-bias off-equilibrium current spin polarization of the homogeneous semiconductor as 

and hence 

P)(x) = P^{x) + P'j{x) ^P + Q^ (136) 

for the total (equilibrium plus off-equilibrium) current spin polarization. 

Turning now to the magnetoresistance of a homogeneous semiconductor, we find from Eq. (|97p that the reduced 
resistance Xh equals that for Ec{x) = A{x) = 0, 

X. - (137) 



[see Eq. (2.69) of Ref. HI]. Therefore, we have, from Eq. p02| . 



R - ^ + 5 

^ h 



for the equilibrium resistance, and, from Eq. (|121p . 



= -2;^K-«2) (139) 



for the off-equilibrium resistance. 
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FIG. 2: Schematic energy diagram for a DMS/NMS/DMS heterostructure composed of three homogeneous layers. 



B. Current spin polarization in DMS/NMS/DMS heterostructures 

For parameterizing a DMS/NMS/DMS heterostructure, we attach labels j — 1,2,3 to quantities referring to the 
left DMS layer, the NMS layer, and the right DMS layer, respectively, and denote the positions of the interfaces 
(including the DMS/contact interfaces) by Xk {k = 1, 2, 3, 4) [see Fig. [2], thereby deviating from the notation used in 
the main body of the paper. Quantities referring to the left and right contact are labelled "f and "r", respectively. 

The position dependence of the polarization inside the left and right contacts (there labelled 1 and 2) is given 
by Eq. (|116p . with fi-{xi^i) = ai^i/p. The current spin polarization Pj{x) in the three layers is found in terms of 
the boundary values of the spin transport function A{x) at the interface positions Xk, using Eq. (|136p together with 
Eqs. (|129p and (|134p . separately in each layer. In each of the layers j, the function A{x) has the Sharvin discontinuities 
/S.Aj and AAj+i at the interface positions Xj and Xj+i, respectively [see Eqs. (3.61) and (3.62) of Ref. jM]. However, 
since in the zero-bias limit, J 0, these discontinuities vanish with J, we may regard the spin transport function A{x) 
as continuous all across the heterostructure. Recalling that in the zero-bias limit the spin-relaxed chemical potential 
/i(a;) may be taken as constant, /i(x) = /ioj it is seen from Eq. (I69p that the function a_(x) is also continuous across 
the entire heterostructure. 

In terms of the quantities ak = a_(xfc) and Ak — A{xk) — e^'^°ak, the complete position dependence of the current 
spin polarization can then be written as 



(0) 



if X < xi and x > X4, respectively, and as 

Pjix)=P,+ 

if Xj < X < I'j+i. Here, 



J. 



(0) 



Gl,r — 



(140) 



(141) 



(142) 



where Gi^r is given by Eq. (|113p . and Qf^ is the Sharvin interface conductance (|10ip . For the coefficients Cji and 
Cj2, we hav( 



,34 



(143) 
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and 



respectively, where 



and 



U{1 - iM, - (1 + 7,)e^^/^^ A,+i] , (144) 



Dj = 2[(1 + 7j) sinh(S'j/Lj) + 27^- cosh(S'j/ij)] (145) 



with Lj and defined by Eqs. (|130p and (|13ip . respectively, and S'j = — a^j is the thickness of layer j. 
Conservation of the total current J across the heterostructure takes the form 

J, = ——jf= J; J = 1,2,3. (147) 

We now require the current spin polarization Pj{x) to be continuous at all interfaces, setting P2 ~ (so that Q2 = 1)- 
This leads to the following set of four coupled linear equations for the quantities ak, 

2 / 

{Gi+Qmhi-Qihia2^ ^Pi, (148) 

VeUl 

2 7 

- Qihiai + {Qigi + e^^''g2) - 6^*^^/1203 = ^ Pi , (149) 

Ven\ 

2 / 

- e^^''h2a2 + (6^*^=32 + Q353)a3 - Qsh^a^ = ^ P3 , (150) 

VeUl 

2J 

- Qshaas + {Gr + Qsga)^^ = 7^ P3 ■ (151) 

Ven\ 



Here, 



and 



gj = hj [cosh{Sj/Lj) + 7j sinh(5j/Lj)] (152) 



h, = ^, (153) 

and we have introduced the band offsets 

S,,, = E^) - . (154) 

The system of equations (|148p - (|15ip is easily solved, so that the complete position dependence of the current spin 
polarization is obtained in explicit form. 

Specializing to the case of a symmetric heterostructure, where the parameters of the right DMS layer and contact are 
identical to those of the left DMS layer and contact, and assuming infinitely high conductivity of the metal contacts, 
tJ/^r 00, so that G/.r — > 00, we find 

ai = 0:4 = , (155) 
2J P ^ ^ 
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and 

2J 

Giai = -Grtti^ 7^P + QhDa2- (157) 

Here, the DMS parameters have been labelled by "D", the NMS parameters by "N", and P = = P^. 

Thus, in the symmetric case, the current spin polarization Pj{x) is completely determined by the quantity 012 • 
Exphcitly, we have, setting L^s'^ — L^P — L^a\ 



Pj{x) = P 



\i X < xi and x > x^, respectively, while 



Pj{x) =P{1- Qg^ + (^^^^hj^ypSoN - xi,i\/LD) + ID sinh(|a; - xi^i\/LD)] \ , (159) 



if xi < X < X2 and X'i < x < Xi, respectively, and 

2PhN 

QgDe~P^'^« +gN + hN 

if X2 < X < X3. Note that, as a consequence of our assumption of infinitely high conductivity of the contacts, the 





Pj{x) - J, f ^— [cosh(^jv/2LAr) + -IN sinh(5Ar/2Lw)] cosh([a; - {x2 + X3)/2]/Ljv) , (160) 

at, as a consequence of our assumption of infinitely high condi 
polarization Pj{x) is independent of the effective electron mass to* and of the densities ri)^ 



C. Magnetoresistance of DMS/NMS/DMS heterostructures 

The total magnetoresistance of a DMS/NMS/DMS heterostructure, R = R + R, is obtained by adding the equilib- 
rium, R, and off-equilibrium, R, contributions corresponding to the three layers [see Eqs. p38p and (|139p ]. 

j=i yj 3 

and 

R=-^^[Piiai~a2) + P3ia3~ai)]. (162) 

For a symmetric structure and infinite conductivity in the contacts, we have 

- 2Q 21d + Sd , 1 21n + Sn i-.nn\ 



and, using Eqs. (|155p and p56p . 



- _ Pa2 _ 1 2P2 



We then have 

R 2P2 



i? + {9N + /iiv)e'3^°« 2Q(2/d + Sd)/21d + e'^^^^ {21n + Sn)/21n 



(165) 



for the relative magnetoresistance R/R = {R — R)/R. As we have assumed Id.n and li^'^'' to be independent of 
the external magnetic field, expression (|165p refiects solely the "spin accumulation part"— of the magnetoresistance. 

''D,N 



While R and R depend on the effective electron mass m* and the densities via the Sharvin interface conductances 
Gi}^N, the relative magnetoresistance R/R is independent of to* and n^^^j^. 
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D. Diffusive limit 

It is instructive to consider in some detail the diffusive limit of the thermoballistic description of spin-polarized 
transport in DMS /NMS /DMS heterostructures as it allows comparison to the results of Ref . 26 obtained within the 
standard drift-diffusion approach to electron transport. 

In the diffusive limit, where Ij < Sj and Ij < we have Lj hi^^ = \fljl^ and 7^ Ij/Li^'' = \/ljJl^ < 1 
[see Eqs. p ^ - ([T^ and pi5|l ]. and hence from Eqs. (fT5^ and (fTOH]) 

gj = hj cosh(S'j/L(^)) (166) 



and 



respectively. Equation (|156p for the symmetric heterostructure then reduces to 



21 ■ 1 

sinh(5,/Li^)) 



"2 = Tm TW\ ' (168) 

Q^Gd coth{SD/Lr') + Gn coth{SN/2Lr') eP^'^^ 



where 



- ^ (169) 



— D^N)^ and Gj is the (spin-summed) conductivity, 

uj = 2(3e^venjlj = TT^f ^J- (I'^O) 

Now, using Eq. (|168p with 6dn — in Eq. (I14ip . we find that the values of the current spin polarization Pj{x) at 
xi, X2, and (2:2 + 2^3)72, respectively, agree with those given by Eqs. (16)-(18) of Ref. [26|, if we identify in the latter 
equations the spin-flip lengths and Aat with l')^^ and Li^\ respectively, and the layer thicknesses d and 2xo with 
Sd and Sn, respectively. 

Similarly, from Eq. (|163p and from Eq. (|164p . using Eq. (|168p . we find the equilibrium contribution to the mag- 
netoresistance, R, and the off-equilibrium contribution, R, respectively, to agree with the corresponding results of 
Ref. viz., the first and second term on the right-hand side of Eq. (10) of that reference. Then, trivially, we find 
agreement of the relative magnetoresistance R/R given by Eq. (|165p with the first term on the right-hand side of 
Eq. (13) of Ref. m. 



E. Ballistic limit 

In the ballistic limit, where > Sj and Ij > we have Lj I'-J'^ and 1 [see Eqs. (|130p - (|132p and (|146p ]. 

so that gj = 1 and hj = e~^^^^'=^ . 

For the symmetric heterostructure, we then have from Eq. (|156p 

2 J P 

^2 — —77^ HTT • (171) 

Qe-/3*°« + 1 + e-SN/e' 

Using this expression for a2, we obtain the ballistic limit for the current spin polarization and the magnetoresistance. 
The equilibrium contribution to the magnetoresistance, in particular, follows from Eq. (|163p in the form 

in which the Sharvin interface resistances l/^/jj^Tv characterize, via their dependence on the constant densities n^^\, 
the (homogeneous) DMS and NMS layers. Notably, owing to the presence of the factor Q, the contribution to R from 
the DMS layers vanishes when P — > 1. If we had not assumed the interfaces a,t x — 2:2,3 to be fixed points of thermal 
equilibrium, the ballistic equilibrium magnetoresistance would be given in terms of the Sharvin interface resistances 
at the end-points xia [see Eq. (2.51) of Ref. |33| and of the potential barriers provided by the profiles ([T]). 
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FIG. 3: Relative zero-bias current spin polarization Pj(x)/P across a symmetric DMS/NMS/DMS heterostructure enclosed 
between metal contacts with infinitely high conductivity, calculated from Eqs. (|158p - (|160|l for zero band offset, 5 = 0. The 
different curves correspond to different values of the momentum relaxation length I = Id = In- The remaining parameter 



values are: Sd = 1 /im, Sn = 2 /im, I 



(N) 



2.5 fim, P = 0.8, L, 



60 nm. 



VI. DMS/NMS/DMS HETEROSTRUCTURES: APPLICATION 

In the thermoballistic description of spin-polarized electron transport in DMS/NMS/DMS heterostructures, the 
momentum and baUistic spin relaxation lengths, I and Ig, respectively, are the fundamental dynamical quantities. 
These, apart from the potential profiles E^^''^^ and the geometric dimensions Sd.N: determine the current spin 
polarization and the magnetoresistance of the structure as a function of the Zeeman splitting A, i.e., of the strength 
of the external magnetic field. Here, we illustrate the dependence of polarization and magnetoresistance on I by 
specific numerical examples. (We recall^i, that, when varying I, we understand the different values to represent a 
class of semiconductors which have similar material properties, but differ in the strength of impurity and phonon 
scattering.) Furthermore, we point out the possibility of determining (experimentally) the quantities I and 1^ from 
their effect on the (experimentally accessible) spin polarization and magnetoresistance. 

In Fig. [21 the relative current spin polarization, i.e., the ratio of current spin polarization and static DMS polariza- 
tion, Pj{x)/P, across a symmetric DMS/NMS/DMS heterostructure at zero bias is shown for different values of the 
momentum relaxation length 1 = Id=In at a, fixed value of the ballistic spin relaxation length Is = li^^ = ls^\ The 
behavior of Pj{x)/P exhibits two features. 

First, the injected spin polarization, i.e., the value of the polarization at the DMS/NMS interfaces, remains on 
a high, weakly Z-dependent level ranging between 0.6 and 0.7 when I varies by three orders of magnitude from the 
diffusive to the ballistic regime. This behavior contrasts with that calculated within the thermoballistic treatment of 
spin injection at the interfaces of the metallic ferromagnets and the semiconductor in a ferromagnet/NMS/ferromagnet 
structure^. The results for the latter case are shown in Fig. 5 of Ref. [13. These have been obtained with an /g-value 
equal to that used in Fig. [3] of the present paper, and with the bulk polarization of the ferromagnets equal to the 
static DMS polarization. The injected polarization at the ferromagnet/semiconductor interfaces (for zero interface 
resistance) is seen to be very small in the diffusive limit, and rapidly rising with increasing I. The strong dependence 
on I reflects the conductivity mismatch^ between metallic ferromagnets and semiconductor, which is very large in the 
diffusive (low-Z) regime where the conductivity of the semiconductor vanishes with I. No such mismatch occurs at the 
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FIG. 4: Relative magnetoresistance R/ R of a symmetric DMS/NMS/DMS heterostructure as a function of the static DMS 
polarization P, calculated from Eq. (|165|l with the parameter values of Fig. [3] (solid curves). The long-dashed curve for / = 10 
nm has been calculated with Sd = 0.2 ^m. 

DMS/NMS interfaces in a DMS/NMS/DMS heterostructure. 

Second, the behavior of the polarization inside the NMS layer is governed by the Z-dependence of the generalized 
spin diffusion length L — [lls/{\ + l/ls)Y^'^ ■ In the diffusive limit, where L w VHI, the polarization dies out rapidly; 
in the ballistic limit, we have L « Zg, i.e., the decay is determined solely by the ballistic spin relaxation length. 
Qualitatively, the behavior of the polarization inside the NMS layer does not depend on the kind of spin injector 
(ferromagnet or DMS). However, due to the larger injected polarization, the magnitude of the polarization inside the 
NMS layer is larger for the DMS injector. 

In Fig. [4] (solid curves), the relative magnetoresistance R/ R of a symmetric DMS/NMS/DMS heterostructure is 
shown as a function of the static DMS polarization P for the parameter values of Fig. [31 For P not too close to 
unity, the qualitative behavior of the curves is determined by the overall factor P^ in Eq. (|165p . When P approaches 
unity, the curves rise sharply owing to the i-*-dependence of the terms proportional to Q in the denominator of the 
expression for R/ R. As a function of Z, R/R rises, over the full P-range, by about one order of magnitude when I 
varies from the diffusive to the ballistic regime. 

In the example of Fig. 21 R/R tends to increase with decreasing thickness Sd of the DMS layers, as is seen by 
comparing the solid and long-dashed curves for Z = 10 nm (the magnitude of this increase becomes progressively 
smaller with increasing Z). In the diffusive calculation of Ref. [2^, the relative magnetoresistance was found, for 
a specific parameter choice, to increase with increasing Sd- We have reproduced this behavior quantitatively by 
choosing the parameter values in Eq. (|165p so as to conform to those underlying the results shown in Fig. 3b of 
Ref. [1^. In particular, we took Z^i = 10 nm and used Eqs. (|132p and (|170p . respectively, to relate the ballistic spin 

relaxation lengths zi^'^^ to the spin diffusion lengths \d,n of Ref. 2&., and the ratio Id/In to the ratio ud/<^n- A 
detailed study of the full parameter dependence of the magnetoresistance is clearly desirable. 

For spin-polarized electron transport in DMS/NMS/DMS heterostructures, the thermoballistic description in the 
simplified form presented in this paper leads to transparent explicit expressions for the current spin polarization 
[Eqs. (|158p - (|160p ] and for the magnetoresistance [Eqs. (|163p and (|164[) ]. If these turn out to reproduce, albeit 
approximately, the results of actual measurements for different values of P and various choices of Sd and 5* at, one 
should be able to determine the values of the parameters Z and Ig for magnetic and nonmagnetic materials from 
experiment in a quite direct manner. 
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VII. SUMMARY AND CONCLUSIONS 



We have presented the systematic extension of the therinobalhstic description of spin-polarized electron transport 
in semiconductors to the case of a spin-split conduction band, allowing us to envisage applications to spin-polarized 
transport in paramagnetic DMS and in DMS/NMS/DMS heterostructures. 

Assuming arbitrarily shaped potential profiles exhibiting arbitrary, position-dependent spin splitting, we have con- 
structed the thermoballistic currents and densities by starting from the spin-resolved densities at points of local 
thermal equilibrium, at which electron currents across ballistic transport intervals are activated. These currents are 
subject to spin relaxation. Dividing the ballistic currents and densities into their equilibrium (spin-relaxed) and 
off-equilibrium parts, wc have obtained a description of the equilibrium parts in terms of the spin-relaxed chemical 
potential, and of the off-equilibrium parts in terms of a spin transport function that is related to the splitting of 
the spin-resolved chemical potentials. From the ballistic currents and densities, the corresponding thermoballistic 
quantities are constructed by weighted summation over all ballistic intervals. The procedures for calculating the 
spin-relaxed chemical potential and the spin transport function are outlined. 

The thermoballistic description of spin-polarized electron transport has been applied to the calculation of the current 
spin polarization and magnetoresistance in DMS/NMS/DMS heterostructures. The results have been compared to 
those of the customary drift-diffusion approach. The important role of the fundamental momentum and spin relaxation 
lengths as well as the possibility of their experimental determination are pointed out. 

The emphasis in the present paper has been on the systematic development of the formalism. In future work, this 
formalism will have to be implemented in detail. To this end, a careful study dealing with the role of the interfaces 
in electron transport in heterostructures, as well as their adequate modelling within the thermoballistic description, 
is required. Furthermore, efficient algorithms for the solution of the integral equations for the resistance function 
and the spin transport function are to be developed. Then, by performing exploratory calculations sampling the 
full parameter space, one may be able to identify novel features and trends in spin-polarized electron transport in 
paramagnetic semiconductors. 
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